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I describe a mechanism to understand the relation between chiral-symmetry breaking and eigenmodes of the 
Dirac operator in lattice QCD with Kogut-Susskind sea quarks. It can be shown that if chiral symmetry is 
spontaneously broken, the eigenvalues Ai should behave as z(i)/V for large volume V, where z{i) is the i-th zero 
of the Bessel function. With neither chiral nor A extrapolation, one can precisely calculate the chiral condensate 
using only a small set of eigenvalues. Therefore, it is economical and free of systematic uncertainties. I present 
the first QCD data to support this mechanism and encourage the lattice community to test and use it. 


One can not completely understand physics of 
hadrons without understanding the QCD vac¬ 
uum, which main properties are confinement and 
chiral-symmetry breaking. 

The chiral order parameter is ('f/’l/') in the chi¬ 
ral limit. Suppose (xx(m, V)) is the one calcu¬ 
lated at bare mass m and lattice volume V. Since 
(xx(0: y)) = 0, conventionally one has to assume 

{tpip) = lim lim {xxijn,V)). (1) 

Practically, most people do extrapolation only on 
finite V and very small set of m by means of some 
unjustified fitting function (linear, quadratic, 
polynomial or logarithm ...). Unfortunately, such 
an process is arbitrary and sometimes gives wrong 
results. Here we provide two known evidences. 

(i) QED in 4 dimensions. Non-compact lattice 
QED experiences second order chiral phase tran¬ 
sition at finite some bare coupling constant g, 
where the chiral condensate vanishes. The critical 
coupling determined by naive m extrapolation of 
(xx(™i^)) is not necessary correct. For detailed 
discussions and current status, see @]. 

(ii) The one-flavor massless Schwinger model. In 
the continuum, as is well known, {ijj'4’)cont = 
0.16e, where e is the electric charge. As indicated 
in Ref. Q , the exact result can not be reached at 
finite bare coupling g = aehy the any extrapola¬ 
tion mentioned above. 

In Ref. [^, we proposed a new mechanism 
to get reliable quantitative information on chiral- 


symmetry breaking, based on the computation of 
the probability distribution function of (iptp)' 

Pic) = lim {6ic - 4 X ^(®)V’(a;)))- (2) 
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Here N is the degrees of freedom, being propor¬ 
tional to V. P(c) is the probability to get the 
value c for the chiral condensate. If there is exact 
chiral symmetry in the ground state, P(c) = S(c). 
If chiral symmetry is spontaneously broken, P(c) 
will be a more complex function. 

For Kogut-Susskind fermions, we derive several 
relations between the chiral condensate and the 
eigenmodes of the Dirac operator. 

(a) The first formula is 
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itptjj) = lim lim ^ / 
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where A^ is the i-th positive eigenvalue of the 
massless Dirac matrix. This is the formula many 
people are currently using. It requires m extrap¬ 
olation, and may not produce reliable result. 

(b) The second formula is 

= lim lim ^(p(A)), (4) 

A —V —»-oo iV 


where p(A) is the eigenmode density. This for¬ 
mula is also known. Since p(0) = 0 at finite V, 
and it requires A extrapolation, which is again 
arbitrary. 



2 


(c) The third formula is 
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( 5 ) 


For its derivation by other methods, see Refs. ^ 
The advantage over the first two formulae is 
that no A extrapolation is necessary. The dis¬ 
advantage is that all the eigenvalues have to be 
calculated. When the lattice volume is large, the 
computational task is huge and not so feasible. 
Indeed, we have checked that it requires a large 
lattice volume for getting consistent results. 

(d) The fourth formula is 


< ipij; >= C{i) 


lim 
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Here z(i) is the i-th zero of the Bessel function 
Jo. This is new and was originally derived in Ref. 
Q. Neither m nor A extrapolation is necessary. 
The advantage over the third formula is that only 
a few smallest eigenvalues are needed for this cal¬ 
culation. The only thing is the hnite size analysis 
which every formula requires to do. 

The relation between eigenmodes and chiral- 
symmetry breaking is clear: if chiral symmetry is 
spontaneously broken, A^ should behave as z{i)lV 
for large V. We have numerically checked Q the 
reliability of Eq. in the Schwinger model, 
where the exact result was reproduced precisely. 
To our knowledge, this is one of the best results 
for the lattice Schwinger model! 

The most interesting application is QCD. Here 
I would like to present the first data for QCD 
with various numbers of quark flavors Nf . In Ref. 
Q, the configurations were generated and eigen¬ 
values of the fermionic matrix were calculated. 
With these, it is very easy to evaluate the chiral 
condensate and compare different methods. 

Figure || shows the quenched result for (3 = 
5.5088 and H = 6^ using the first formula: Eq. 
(^). Only a small set of data can be fitted by 
a linear function. If one hts only the data for 
m £ [0.04,0.08] to the chiral limit, the result is 
0.060. But outside this range, either for larger or 
smaller m, this slope changes and global fitting 
function can not be easily found. Eor smaller m, 


the change is very rapid. We attribute it to finite 
size effects. 

Figure || shows the quenched result for [3 = 
5.5088 and V = 6"*, but using the fourth formula: 
Eq. (H). Even on such a small lattice, there is a 
nice plateau for i £ [10,25], from which one can 
reliably get C{i) = 0.06. No fitting function is 
necessary. Of course, detailed finite size analysis 
remains to be done. 

Figures I I i I are the results for Nf = 
1, 2, 3,4. All these satisfy well Eq. (^). 

In conclusion, I have shown a reliable method 
for investigating the chiral properties and obtain¬ 
ing the chiral condensate from only a small set of 
A, without m or A extrapolation. This method 
was previously tested in the Schwinger model. 
In this paper I have presented the first data for 
QCD. Therefore, I encourage other people in the 
lattice community to test and use Eq. (^. 
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Figure 1. ^)) for Nj = 0, f3 = 5.5088 and 

V = 6'^ using Eq. (3). 
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Figure 2. C{i) as a function of i for Nf = 0, 
(3 = 5.5088 and E = 6^ using Eq. (6). 
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Figure 3. C{i) as a function of i for Nf = 1, 
P = 5.3469 and E = 6^ using Eq. (6). 
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Figure 4. C{i) as a function of i for Nf = 2, 
P = 5.3136 and E = 6'^ using Eq. (6). 


0.1 
0.095 
0.09 
0.085 
C{i) 0.08 
0.075 
0.07 
0.065 
0.06 

10 15 20 25 30 35 40 

i 

Figure 5. C{i) as a function of i for Nf = 3, 
P = 5.2814 and E = 6^ using Eq. (6). 
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Figure 6. C{i) as a function of i for TV/ = 4, 
(3 = 5.2500 and V = using Eq. (6). 




